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摘 要 : 首先 给 出 了 不 可 约 非 负 人 矩阵 最 大 特征 值 的 上 上 下界。 然后 利用 相似 变换 构造 了 一 列 相似 矩阵 ， 从 而 
得 到 不 可 约 非 负 和 矩阵 最 大 特征 值 的 逐步 压缩 的 一 列 上 下 界 ， 其 极限 为 所 要 求 的 最 大 特征 值 。 最 后 
利用 2G- 矩 阵 与 非 负 和 矩阵 的 关系 ， 给 出 了 计算 不 可 约 2- 和 矩阵 最 小 特征 值 的 一 个 新 算法 。 理 论 上 给 
出 了 收敛 性 证 明 。 该 算法 迭代 过 程 简 单 ， 不 用 计算 逆 矩 阵 ， 从 而 计算 量 小 ， 占 用 内 存 少 。 数 值 实 
验 的 结果 表明 该 算法 具有 可 行 性 和 有 效 性 。 
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Ze" = {A = (aij) € R™” | aij <0, 1<i, j <n, if#j}. 


进一步 ,对 B = (b;)>0, A= sI-B, BAEZ”, Hs > p(B) WRAAA 
异 M- 和 矩阵 。 若 对 s CIBER, B = (bij) > 0， 则 称 矩 阵 4 为 2- 矩阵。 显然 ，2G- 矩 阵 是 
比 M- 和 矩阵 更 广泛 的 矩阵 ， 出 现在 许多 应 用 领域 中 ， 经 常 需要 计算 2- 和 矩阵 的 最 小 特征 值 ， 而 通 
常 的 方法 是 先 对 2- 矩阵 的 上 下 界 进行 估计 ， 见 文献 [1-11]。 文 献 中 基 于 他 们 给 出 的 2- 矩阵 最 
小 特征 值 上 下 界 的 一 个 估计 ， 提 出 了 求解 3F- 矩 阵 最 小 特征 值 及 特征 向 量 的 数值 算法 ， 但 他 们 方 
法 的 收敛 速度 较 慢 。 本 文 利用 2- 矩 阵 与 非 负 矩阵 的 关系 ， 给 出 了 计算 不 可 约 2 矩阵 最 小 特征 
值 的 一 个 新 算法 。 该 算法 迭代 过 程 简单 ， 与 文献 [7] 和 文献 [8] 相 比 ， 收 敛 速度 更 快 。 
设 G= (gij nxn 是 不 可 约 2 和 矩阵， 令 


R > max |gs], tEN, N = {1,2, ,n}. 
于 是 有 A = RI 一 G = (aij)axn HP Ln POs, WN A AAA EXN AJ (an > 0) 的 非 
负 不 可 约 矩 阵 。 如 果 Xo 是 4 的 最 大 特征 值 ， 则 w = R 一 和 0 为 G 的 最 小 特征 值 。 因 此 ， 不 可 
约 2- 矩 阵 G 的 最 小 特征 值 的 算法 依赖 于 具有 正 对 角 元 的 非 负 不 可 约 矩 阵 44 的 最 大 特征 值 的 算 


2 ”主要 结果 
首先 给 出 以 下 两 个 引 理 ， 参 见 文献 [6]。 
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引 理 1 Bai, ap, ary bi,be,---, b 是 任意 的 实数 ， 则 


i ai ai +02 +- +an i Qj 
当 且 仅 当 所 有 的 
t= 1,2,- ‘Nn 
ay 


相等 时 ，(1) 式 等 号 成 立 。 
引 理 2 设 4= (aij)nxn，A* 的 第 i 行 行 和 与 第 i 列 列 和 分 别 记 为 ri(A*) 与 ci(A*)， 则 


ri(A*+!1) = ary ci(A*+1) = Sat (2) 
t=1 


以 下 总 记 AW n ARTA AES. 
定理 1 设 入 是 4 的 最 大 特征 值 ，B = (4 十 了 )? (7 为 单位 矩阵 ，p 为 任意 给 定 的 正 整 数 ) 。 


以 mi(C4) 和 ci(4) 分 别 记 妇 的 第 ; 行 行 和 与 第 ; 列 列 和 。 若 mi(4) > 0, i 二 1… ,n， 则 


— < < 
me (B) < ìo S max r:(B) ’ (3) 


#al(A)>0, i=1,…,n， 则 


.Ci(AB) Ci 
< < ; 
ey SD 


证 明 类 似 于 文献 [6] 中 定理 1 的 证 明 。 
记 Ao =A= (Qij jnxn? 构造 如 下 和 迭代 序列 
Ák+1 = Dy Ax Dy, k=0,1,2,---, 
其 中 
Dk = diag(rı (Ax), r2(Ax), we ,Tn(Ak)); 


显然 矩阵 序列 {An} 是 相似 的 且 不 可 约 。 又 记 Bk = (A, +1)? (p 为 任意 给 定 的 正 整数 )。 于 是 可 
得 本 文 的 主要 结果 。 


R (k) in ikB) p) ri( Ak Br) 
Tmin = CRB Tmax = MAX — Be)” 
WW A 的 最 大 特征 值 和 o 满足 
AO arD eiar er (5) 
证 明 记 


(A,B 
(6 men = Ansa Aon +P, (dG wan = (Aner EP, fO = BEE) 


则 


+1 rj(Ak) (k E k 
dD Eud È rjad 


ze l 
: Sa =: Aaa} D ri(An)ayy” 
j= j= j= 
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由 引 理 1 有 

n Ck) ~ (Kk) 

(k) È> a pet) 名 ig 

Tmin 一 TE min n w Í 一 Ti < max n (k) = re. (6) 
2 dij X dij 
j=l j=l 
即 rr) < hth) < riy 从 而 

eS ee 


(7) 
下 面 的 定理 保证 了 区 间 序 列 rE, re 


min? Tmax |; k = 0, 1,: ee 将 逐步 收 罕 。 
定理 3 ”对 任意 迭代 步 上， 存在 一 个 正 整数 s (1 < s < n)， 使 得 


(kts) < rE), rikt+s) > re. 


证 明 类 似 于 文献 7] 中 定理 2 的 证 明 。 
与 文献 7] 比较， 本 文 有 关 不 可 约 非 负 和 矩阵 的 最 大 特征 值 的 算法 同样 是 用 对 角 变 换 进行 办 
代 ， 在 迭代 过 程 中 对 角 元 始终 保持 不 变 ， 非 零 非 对 角 元 的 计算 只 需 采 用 简单 的 公式 c(*+1) = 


MOO /Jr ， 因 此 计算 量 不 大 ， 占 用 内 存 少 ， ASE SOAS HO 在 保 
持 文 献 中 算法 优点 的 同时 ， 本 文 算 法 的 收敛 速度 更 快 。 我 们 有 如 下 结 
定理 4 对 任意 的 正 整 数 m (m = 0,1，: 


oe 
r(A(A+1)™*1) ri(A(A+D)™) 
(A+ Dm) > PX 


i ri((4 + 1)®) i 
_ ri(A(A + 1)™*1) . Ti(A(A+7)™) 
mn aD) > wh AT 

证 明 WAHI = (bij)nxn 由 引 理 1 和 引 理 2 得 


n 


ri(A(A + 1)™+1) 当 | para) 
n(A +N) ~ 


> biere(A + I)™ 


Sg ae a 
Se 
min ACA +) So BAA +D") 


ARO = cae = 
当 我 们 给 定 适当 的 正 整 数 m 时， 每 进行 一 次 对 角 变 换 ， 利 用 本 文 方法 估计 得 到 的 最 大 特征 值 上 


下 界 ， 其 界 宽 比 文献 四 方法 的 界 宽 窗 。 因 此 ， 本 文 方法 的 收敛 速度 明显 比 文献 7 的 方法 快 
达到 给 定 精度 所 需要 的 迭代 次 数 更 少 。 


i (At Dm) =e 
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于 是 ， 基 于 不 可 约 非 负 矩阵 与 不 可 约 2 矩阵 的 关系 ， 我 们 有 下 面 的 结果 。 
定理 5 设 G 是 不 可 约 ZAR, w ÆG KRIM IEH. MR>O0, 使 RI 一 G = A 和 4 是 不 可 约 
非 负 失 阵 ，4 的 最 大 特征 值 Xx， 则 ww = R- 和 Xo。 


3 ”算法 及 数值 实验 

算法 

步骤 0 输入 不 可 约 2- 和 矩阵 G = (gij)nxn, €> 03 

步骤 1 取 R=1+max|gil， S Ao = RI -G = (aij)nxn, Bo = (ho 十 I)?， 了 为 n 阶 单位 
和 矩阵，p 为 任意 给 定 的 正 整数 ，k = 0; 

步骤 2 计算 


(k) _ Ti(ARBr) (k) 
Ti we ri(Bx) 》 Tmin 


(k). 


i 9) 


(k) 


(k) — 
; Tmax = maxr 


= minr. 
i 


步骤 3 WEE- rE) < es， 则 转 步 骤 5; 
步骤 4 +> 


D; = diag(r1(Ax),r2(Ak),*** ,rn(Axr)), Ak+1 = DE AkpDE, 
Bk+1 = (Angi +I», k—k+1, 
LEAD RR 2; 


步骤 5 MAw=R- Hre +rE) 
例 1 设 


G 为 不 可 约 Z- K G 的 最 小 特征 值 。 
例 2 设 
8 0 -2 -1 0 -1 0 0 


G 为 不 可 约 2 矩阵 ， 求 G 的 最 小 特征 值 。 
利用 本 文 给 出 的 算法 ， 我 们 可 以 求 得 例 1 和 例 2 算 阵 的 最 小 特征 值 ， 并 与 文 了 7] 和 文 [8] 的 算 
法 比较 ， 数 值 结果 分 别 见 表 1 和 表 2。 
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表 1: 例 1 的 计算 结果 及 比较 


[7] 的 选 代 次 数 本 文 算法 PP = 三 5 时 的 迭代 次 数 精度 w 
13 5 1074 -4.7736 
25 16 1078 -4.77363511 
31 22 10710 -4.7736351183 
36 28 107}? -4.773635118400 


表 2: 例 2 的 计算 结果 及 比较 


[7,8] 的 迭代 次 数 。” 本文 算法 p= 二 2 时 的 迭代 次 数 HE w 
17 15 1074 0.9444 
33 30 1078 0.94440470 
48 46 107+? 0.9444046950296 
67 64 10776 0.9444046950294975 
67 64 10-20 0.94440469502949753178 
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A Numerical Algorithm for the Minimal Eigenvalue of 
an Irreducible Z-matrix 
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Abstract: It is well-known that the problem to calculate the minimal eigenvalue of a matrix commonly 
occurs in many branches of science and engineering. In this paper, we consider the problem to calculate 
the minimal eigenvalue of a so-called Z-matrix. An upper and a lower bounds on the maximal eigenvalue 
of a nonnegative matrix are firstly given. Then by using a similarity transformation, a series of upper 
and lower bounds on the maximal eigenvalue of a nonnegative matrix are obtained, these bounds 
gradually approach the maximal eigenvalue. Finally, based on the relation between the Z-matrix and 
the nonnegative matrix, a numerical algorithm for computing the minimal eigenvalue of an irreducible 
Z-matrix is proposed. It is shown by numerical examples that the convergency rate of the presented 
method is faster than that of previous methods. 

Keywords: nonnegative matrices; Z-matrices; irreducible; minimal eigenvalue; convergency rate 
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